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Abstract:  The  intent  of  this  research  is  to  study  the  dynamic  behavior  of  a 

solid  body  resting  on  a  moving  surface.  Results  of  the  study  are  then  used  to  pro¬ 
pose  methods  for  controlling  the  orientation  of  parts  in  preparation  for  automatic 
assembly.  Two  dynamic  models  are  discussed  in  detail.  The  first  examines  the 
impacts  required  to  cause  reorientation  of  a  part.  The  second  investigates  the  use 
of  oscillatory  motion  to  selectively  reorient  parts.  This  study  demonstrates  that 
the  dynamic  behaviors  of  solid  bodies,  under  the  conditions  mentioned  above,  vary 
considerably  with  small  changes  in  geometry  or  orientation.  > 
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1  Introduction  and  Problem  Definition 

-‘(Conventionally,  parts  are  oriented  by  bowl  feeders.  These  machines  vibrate  and 
thereby  convey  parts  through  a  series  of  filters  which  reject  all  but  a  particular  orientation. 
Rejected  parts  are  returned  to  storage.  Many  researchers  have  examined  the  implemen¬ 
tation  of  programmable  or  adjustable  filter  stages^  (Boothroyd,  1975)  (Boothroyd  and  Ho, 
1977)  (Boothroyd  and  Vlurch.  1970)  (Boothroyd  et  ah,  1982)  (Boothroyd  et  al.  1977) 
(Lozano-Perez,  1986)(Murch,  1977)  (Murch  and  Boothroyd,  1975)  (Murch  and  Poli,  1977) 
(Redford  et  al.,  1983a)  (Redford  et  al.,  1983b)  (Singer,  1985)  (Suzuki  and  Kohno,  1981). 
These  techniques  are  often  successful  for  limited  classes  of  parts.  The  scope  of  this  paper 
is  to  present  theoretical  results  which  may  be  useful  in  feeder  designs.  It  is  hoped  that  a 
more  detailed  understanding  of  t  he  dynamic  behavior  of  bodies  resting  on  a  moving  surface 
will  facilitate  the  design  of  new  types  of  programmable  parts  feeders. 

^  ■'  <•.  /'  '  -  ‘  If  /'■  «'c  ftc  ^  \ 

•  - 

2  Impact  Reorientation 

In  this  section,  we  determine  the  conditions  necessary  to  cause  part  reorientation.  A 
part  starts  in  a  stable  orientation  (a  natural  resting  aspect  (Boothroyd  and  Ho.  1977))  on  a 
flat  surface.  It  is  then  given  an  initial  horizontal  velocity  after  which  it  impacts  a  wall  near 
the  surface  upon  which  it  is  resting  (figure  1).  The  equations  describing  the  impact  are 
based  on  conservation  of  angular  momentum  about  the  impact  point.  After  impact,  the 
principle  of  conservation  of  energy  is  applied  in  order  to  determine  whether  reorientation 
occurs. 

The  part  shown  in  figure  1  is  a  cylinder,  but  the  analysis  easily  generalizes  to  any 
object.  The  important  parameters  are  the  mass,  m  of  the  part;  the  moment  of  inertia, 
I :  and  the  constants,  /  and  h.  which  locate  the  center  of  mass,  cm,  with  respect  to  the 
impact  point.  The  expressions  for  angular  momentum,  L,,  before  and,  L/,  after  impact 
are  (Klepner  and  Kolenkow,  1973): 

h i  -  mVlh 

,  __  IgVf  (2-1) 

'f  Vh-  ■  l2 

“hi. re  '  is  the  speed  of  the  center  of  mass  prior  to  impact,  V /  is  its  speed  just  after  impact, 
and  the  inertia  about  the  impact  point,  is  defined  as 
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the  change  in  potential  energy  which  will  cause  the  part  to  reorient  yields 
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Solving  equations  2.1  for  Vf,  subst  itut  ing  into  equation  2.2,  and  t  hen  solving  the  result 
for  the  critical  value  of  F,  yields 

-  I'1)'  *1*  ,(*’  '  I1)  ■ 

V„a  -  - 1 - 1 -  .  (2.15) 

which  is  an  equation  for  the  minimum  initial  horizontal  velocity  that  will  cause  reorienta¬ 
tion.  For  cylinders,  ^  is  ^ h 1  +  |/2;  for  rectangular  parts,  is  z{h'i  f  l2). 

In  order  to  test  the  ability  of  impact  reorientation  to  discriminate  among  different 
orientations,  it  is  useful  to  consider  a  part  which  has  two  very  similar  orientations.  Figure 
2a  shows  a  rectangular  solid  part  in  two  possible  orientations.  One  edge  is  bb  longer 
than  the  other  where  ~  is  expressed  as  a  fraction  less  than  one.  As  the  part  exchanges 
orientations,  the  center  of  mass  location  changes  in  both  height  and  horizontal  distance 
from  the  impact  point  by  6.  The  ratio  of  F,ri(  of  orientations  one  and  two  (denoted  by  F( 
and  V'2)  expanded  in  a  power  series  in  6  is  expressed  as: 

F, 

~  -  I  +  2.2 b  •  1.336“  -  •  .  (2,1) 

n 

This  equation  illustrates  that  small  differences  in  location  of  the  center  of  mass  result  in 
significant  variations  in  impact  reorientation  velocity. 

To  verify  the  equations  and  to  get  a  measure  of  experimental  error,  tests  were  per¬ 
formed.  Cylinders  of  different  size  were  placed  on  a  moving  conveyor  belt.  The  cylinders 
impacted  a  steel  wall  .889  mm  (.035  inch)  high.  Fifteen  trials  were  made  for  each  velocity 
set  point  and  the  number  of  pegs  which  reoriented  was  counted.  The  velocity  of  the  belt 
was  monitored  by  a  tachometer  which  was  friction  driven  directly  by  the  belt.. 

Several  tests  were  performed.  First,  the  lower  threshold  velocity,  Fy^,  was  determined 
by  finding  the  maximum  velocity  at  which  no  pegs  fell  over  in  fifteen  trials.  The  upper 
threshold  velocity,  \’jaih  was  determined  similarly  by  finding  the  lowest  velocity  at  which 
all  pegs  fell  over.  The  median  velocity  Vme,i  is  simply  the  center  of  the  velocity  band 
defined  by  Vj^ji  and  F/„//.  Finally,  two  percentages  have  been  determined  from  this  data. 
The  percentage  of  error  is  calculated  between  the  theoretical  velocity,  VT^rory,  and  the 
median  experimental  velocity,  Vnieil.  The  uncertainty  band  of  the  experiment,  is  calculated 
from  and  Vf„u\  this  band  is  represented  as  a  percentage  of  Fm,.,p 

Results  confirm  that  the  system  model  adequately  represents  the  experimental  situ¬ 
ation.  All  tests  performed  had  narrow  uncertainty  bands.  This  indicates  that  there  is  a 
sharp  cutoff  between  not  having  enough  energy  to  cause  reorientation  and  having  enough 
energy  to  guarantee  reorientation.  In  addition,  the  error  percentages  between  Fr/,.  ,,,,  and 
F,,,, are  small.  Fur* hormoro,  the  scaling  ol  the  velocity  agrees  with  the  scaling  predicted 
in  equation  2.3  Two  pegs  which  vary  greatly  in  size  and  length  to  diameter  ratio  can  have 
the  same  reorientation  "elocity.  Table  1  gives  typical  results  for  two  different  cylinders. 

The  ability  to  discriminate  between  orientations,  or  to  reorient  one  object  while  not 
reorienting  another,  depends  heavily  on  the  quality  of  the  velocity  source.  Fortunately, 


excellent  velocity  sources  are  inexpensive  and  easy  to  build.  Figure  2b  presents  a  velocity 
line  with  the  reorientation  velocities  for  the  two  parts  shown  in  figure  2a.  The  dark  section 
of  the  line  identifies  the  desired  operating  velocity  range.  Within  this  range  part  1  will 
always  be  reoriented,  while  part  2  will  never  be  reoriented.  By  cascading  several  of  these 
impacts  at  different,  velocities,  a  set  of  like  parts,  initially  oriented  randomly  can  all  be 
driven  to  the  same  orientation. 


3  Vibratory  Reorientation 
3.1  Modeling 

The  model  for  this  section  is  an  object  in  a  gravity  field.  The  table  on  which  the 
object  rests  oscillates  vertically  (Figure  3).  The  object  is  first  given  some  initial  angl  e,  On, 
possibly  by  a  sudden  motion  of  the  table.  For  a  proper  choice  of  X(t),  the  table  driving 
function,  the  given  object  can  be  maintained  in  a  rocking  mode  of  fixed  amplitude.  The 
driving  function  of  the  table  can  then  be  changed,  causing  an  increase  in  the  amplitude  at 
which  the  object  will  continue  rocking.  Finally,  the  amplitude  can  be  increased  until  the 
object  reorients. 

Because  of  the  discontinuity  of  motion  as  the  part  impacts  the  table,  the  solution  to 
the  equations  governing  the  system  must  be  broken  into  regions.  First,  the  region  in  which 
there  is  no  impact  will  be  examined.  The  equations  of  this  motion  may  be  generated  using 
Lagrange's  method.  The  potential  energy  function  can  be  expressed  as 


S 
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V  ~  mg(r  sin  0  +  x)  ,  (3.1) 

where  V  is  the  potential  energy,  m  is  the  part  mass,  r  is  the  distance  from  the  contact 
point  to  the  center  of  mass,  0  is  the  angle  between  the  horizontal  and  the  vector  r,  and  x 
is  the  vertical  position  of  the  table.  The  kinetic  coenergy  can  be  expressed  as: 


T-  lr  I  d<t>\2  1  \dr  dd>  1  f  d<t>  ■  /o  9\ 

2  \dt  )  2  [dt  dt  J  2  [  dt  J  v  ’ 

where  lc  is  the  moment  of  inertia  of  the  part  about  the  center  of  mass.  These  two  terms 
can  be  used  to  form  the  Lagrangian  which  is  then  substituted  into  Lagrange’s  equation, 
yielding  the  two  equations  of  motion  for  the  system, 
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dt'  \ 


"i-—  -  mr-—  cos O  mr—  sin  0  F  mg  —  F  , 
dt'  dt*  dt 

The  first  of  the  equations  gives  the  basic  motion  of  the  object,  while  the  second  equation 
is  an  expression  for  F,  the  contact  force.  These  equations  are  valid  in  the  conservative 
regions  without  impact. 


The  next  step  in  formulating  the  model  is  to  represent  the  impact  of  a  part  with  the 
table.  At  impact,  the  center  of  mass  of  the  object  has  a  vertical  velocity  which  is  reflected 
with  a  coefficient  of  restitution,  e.  In  a  perfectly  elastic  impact,  f  is  equal  to  one.  Physical 
experiments  were  performed  to  determine  the  coefficient  of  restitution  for  aluminum  parts 
rocking  on  an  aluminum  plate.  The  coefficient  of  restitution  was  determined  to  be  between 
.7  and  .75  for  these  tests.  For  all  of  the  simulations  in  this  paper,  an  e  of  .5  was  used  to 
be  conservative. 

As  a  part  hits  the  table,  the  impact  will  also  affect  the  horizontal  velocity.  If  the  table 
is  frictionless,  the  part  will  slide.  In  this  derivation  the  table  was  assumed  to  be  rough 
enough  so  that  the  contact  surface  will  not  slide  on  the  table  and  after  impact  the  part 
w'ill  continue  to  rock  by  rotating  about  the  new  contact  point. 

The  next  step  in  creating  steady  rocking  motion  was  to  determine  table  acceleration 
profiles  which  would  cause  a  part  to  oscillate  at  fixed  amplitude.  The  amplitude  of  vibration 
is  coupled  with  the  frequency  of  oscillation.  As  the  part  gains  energy,  both  the  amplitude 
and  the  period  of  oscillation  become  greater.  This  response  makes  analysis  difficult:  few 
tools  are  available  for  predicting  the  behavior  of  this  type  of  system  (Meirovitch,  1975). 

An  analytical  solution  to  the  system  equations  for  the  time  required  for  passage  from 
an  initial  to  a  final  angle  can  be  derived  for  the  case  in  which  the  driving  function,  ^4 
remains  constant.  The  derivation  of  this  solution  is  based  on  conservation  of  energy  in 
the  region  without  impact.  For  a  constant  table  acceleration,  the  part  can  be  considered 
to  be  in  a  conservative  field  as  shown  in  equation  3.3a.  In  this  particular  case  the  energy 
balance  equation  becomes 
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where  d>n  is  the  starting  angle  of  the  object,  and  <j>  is  some  angle  of  interest.  Note  that  the 
derivative  of  this  expression  is  the  result  obtained  in  equation  3.3a.  Equation  3.1  can  be 
solved  for 
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This  equation  can  in  turn  be  solved  by  separation  of  variables  and  then  once  integrat¬ 
ing  with  respect  to  time.  For  convenience  the  constant  A  is  defined  as 
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and  the  result  is 
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where*  o..  is  the  ending  angle  of  interest .  The  right  side  is  a  simple  integral;  the  left  side  is  an 
elliptic  integral.  The  solution  to  this  equation  will  be  presented  as  a  combinat  ion  of  elliptic 
functions.  Though  elliptic  functions  can  not  be  expressed  in  terms  of  elementary  functions, 
they  are  well  understood  and  are  easih  computed.  In  order  to  show  that  equation  11.7 
becomes  an  elliptic  integral,  a  change  of  variables  must  be  made,  bet 
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(3.H) 


Substituting  equations  3.11  and  3.8  into  the  left  hand  side  of  equation  3.7  produces  the 
following  expression 
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In  the  case  of  impact 
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where  ij>  is  defined  as  the  angle  between  the  bottom  of  the  part  and  the  vector  r  from 
the  contact  point  to  the  center  of  mass.  If  the  integral  is  then  rewritten  as  two  integrals 
starting  from  zero,  equation  3.12  can  be  expressed  in  terms  of  the  elliptic  integral  of  the 
first  kind  (Gradshteyn  and  Ryzhik,  1080): 
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The  elliptic  integral  of  the  first,  kind  expressed  in  normal  trigonometric  form  is 


(3-13) 
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in  terms  of  the  system  parameters. 

3.2  Choice  of  Driving  Function 

As  mentioned  previously,  this  analysis  is  correct  only  for  constant  values  of  (^4.  The 
first  driving  function  to  be  considered  was  a  periodic  rectangular  wave  in  acceleration  (fig¬ 
ure  4).  Four  parameters  must  be  selected  to  fully  specify  the  waveform  the  acceleration 
upward,  AXl;  the  acceleration  downward.  AX2:  and  the  switch  times.  Tl  and  T2.  The 
selection  of  these  parameters  is  coast  mined  by  several  relations.  First,  the  areas  AJ  and 
A2  must  be  equal  in  order  for  the  table  to  return  to  its  initial  position  and  velocity  after 
each  cycle.  Second,  the  part  must  start  and  finish  the  cycle  at  the  same  angle  (the  sta¬ 
bility  requirement).  Lastly,  the  acceleration  downward  must  never  exceed  gravity.  This 
constraint  is  set  so  that  the  part  does  not  lose  contact,  with  the  table. 

Equations  3.13  and  3.15  were  substituted  into  a  nonlinear  solution  program  which 
calculated  accelerations  and  switch  times  that  brought,  a  chosen  part,  back  to  its  initial 
angle  in  exactly  one  cycle.  Figure  5  shows  a  time  history  of  the  acceleration  of  the  table, 
jjjj-,  and  of  the  angle  (j>  as  determined  by  numerical  simulation  of  equation  3.3  including  the 
model  for  impact.  For  this  case,  the  rocking  motion  is  stable.  However,  small  parameter 
variations  cause  an  instability.  Figure  6  shows  the  rocking  motion  for  an  object  given  a 
slightly  larger  initial  amplitude  than  that  of  figure  5.  After  the  first  impact,  the  peak 
amplitude  is  too  low;  after  the  second,  it  is  too  high.  Eventually,  the  phasing  between 
input  and  output  is  lost  and  the  part  ceases  to  rock. 

Because  our  goal  is  to  bring  parts  to  a  stable  rocking  motion  from  a  range  of  ini¬ 
tial  conditions,  the  rectangular  acceleration  waveform  will  not  do  the  job.  A  variety  of 
waveforms  were  tested.  One  which  exhibited  good  performance  is  shown  in  figure  7.  The 
waveform  is  simply  the  rectangular  waveform  with  the  edges  sloped.  Figure  8  shows  results 
from  a  numerical  simulation  using  this  forcing  function  as  input.  A  part  is  started  with 
an  initial  angle  which  is  twice  the  steady  state  angle  yet  it  still  reaches  a  stable  equilib¬ 
rium  motion.  This  suggests  that  by  using  the  proper  waveform,  a  great  range  of  initial 
conditions  of  the  part  can  be  tolerated. 
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If  this  method  is  to  bo  used  to  selec  tively  reorient  parts,  different  parts  should  exhibit 
significant ly  different  responses.  Figure  9  shows  variation  in  response  for  two  similar  parts 
given  the  same  initial  angle  and  subjected  to  the  same  acceleration  profile.  For  these 
simulations,  upward  acceleration  was  set  to  occur  around  the  time  of  impact.  Several 
simulations  were'  executed  with  various  length  "windows”  of  upward  acceleration.  The 
ability  to  discriminate  degraded  as  the  window  became  very  large  in  comparison  with  the 
difference  in  impact  times. 

An  attempt  was  made  to  design  an  experimental  apparatus  for  study  of  the  behavior 
described  above.  No  satisfactory  configuration  could  be  created  at  a  reasonable  cost.  To 
drive  reasonably  large  parts  with  periods  of  oscillat  ion  greater  than  .5  seconds  would  require 
a  table  with  a  range  of  motion  substantially  greater  than  that  of  available  electromagnetic 
shakers.  Hydraulic  shakers  are  designed  to  deliver  large  forces  and  have  servo  valves  which 
are  too  small  to  supply  the  flow  rates  necessary  to  achieve  the  required  velocities.  The 
ideal  driving  system  would  be  a  small  diameter  hydraulic  actuator  with  a  relatively  large- 
servo  valve. 

4  Conclusion 

This  study  of  part  motion  has  demonstrated  that  small  variations  in  inertial  propert  ies 
from  one  part  to  another  cause  significant  differences  in  the  dynamic  behavior  of  the  two 
parts.  Two  example  techniques  which  capitalize  on  this  property  were  presented.  Practical 
implementation  of  the  first  technique,  impact  reorientation,  is  feasible  given  current  design 
technology.  The  second  technique,  vibratory  reorientation,  will  require  construction  of 
oscillating  tables  with  performance  characteristics  available  only  through  custom  design. 
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Fig.  5a.  Response  of  a  part  to  input  in  fig  5b. 


